Method for blind channel estimation

ABSTRACT

Provided is a method for channel estimation increasing frequency band efficiency lost by using a pilot, and reducing sensitivity to channel zero, instability, and complexity. The method includes the steps of: (a) generating an i-th symbol block S i  including N carriers; (b) performing an inverse fast Fourier transform (IFFT) operation on the i-th symbol block and forming an orthogonal frequency division multiplexing (OFDM) symbol block; (c) attaching a guard interval sample in front of the i-th OFDM symbol block U i  and forming at least one OFDM symbol block U i,cp ; and (d) modeling the formed OFDM symbol block U i,cp  with a channel finite impulse response (FIR) filter h and noise v and estimating channel impulse response using signals yi received through a channel.

CROSS-REFERENCE TO RELATED APPLICATION

This application claims priority to and the benefit of Korean Patent Application Nos. 2005-119519, filed Dec. 8, 2005, and 2006-69472, filed Jul. 25, 2006, the disclosures of which are incorporated herein by reference in their entirety.

BACKGROUND

1. Field of the Invention

The present invention relates to an orthogonal frequency division multiplexing (OFDM) system, and more particularly, to a method for blind channel estimation using a guard interval of a signal in the OFDM technique.

2. Discussion of Related Art

In order to provide telecommunication services in various forms such as voice and packet data, wireless communication systems are widely disposed. The systems may be multiple access systems that share available system resources to support communication with multiple users. A code division multiple access (CDMA) system, a time division multiple access (TDMA) system, and an orthogonal frequency division multiple access (OFDMA) system are examples of multiple access systems.

Here, the OFDM is a technique that separates data into several bit streams and modulates the bit streams using several carriers. More specifically, the OFDM technique converts a serial bit stream into a parallel bit stream and modulates and transmits the information using subcarriers of different frequencies. In comparison with a method transmitting data in sequence using one carrier, the OFDM technique has a longer interval between transmission signals and thus is less affected by a channel delay time. Also, the OFDM technique can reduce interference between successive signals and thus is strong for a multipath channel. In addition, the OFDM technique can increase spectrum efficiency and shows good bandwidth efficiency because high-speed transmission is possible.

Therefore, the OFDM technique is hardly affected by time delay of a multipath and thus does not require a time-domain equalizer. In addition, the OFDM technique can eliminate interference between signals by inserting guard intervals.

In general, according to the OFDM technique, when there is no adjacent signal interference or adjacent channel interference caused by distortion of a transmission channel, orthogonality between subchannels is maintained, and each subchannel is completely separated by a fast Fourier transform (FFT) operation at a receiving end.

However, since a spectrum of an OFDM signal is not a limited band in actuality, energy of each subchannel is transferred to an adjacent channel due to linearity distortion such as multipath. Thus, adjacent signal interference is caused. While this problem can in theory be solved by increasing the number of carriers or a signal period, in practice this is difficult due to carrier stability, Doppler shift, and FFT size. Thus, instead, guard intervals are inserted into an OFDM signal.

FIG. 1 is a diagram illustrating a method for estimating a blind channel by inserting guard intervals into an OFDM system according to conventional art. Each OFDM symbol consists of two parts, i.e., an actual signal interval 10 and a guard interval 20. Here, the guard interval 20 is made by attaching an end part of a signal to a beginning part thereof.

Referring to FIG. 1, assuming that a k-th carrier in an i-th symbol block is s_(i)(k) in OFDM signal transmission, an i-th symbol block consisting of N carriers is expressed by a column vector s_(i) =[s_(i)(0), . . . , s_(i)(N−1)]^(T).

Here, each symbol s_(i)(k) is statistically uncorrelated with the others as in most known modulation methods (binary phase shift keying (BPSK), phase shift keying (PSK), quadrature amplitude modulation (QAM), differential phase shift keying (DPSK), and so on), independent of the others, identically distributed, and has a zero average and a variance of 1.

Here, in the OFDM system, the i-th symbol block s_(i) is modulated into a time-domain signal 10 by an N-point inverse fast Fourier transform (IFFT) operation. The time-domain signal 10 is expressed by the following Formula 1. u _(i) =[u _(i)(0), . . . , u _(i)(N−1)]^(T) =F _(N) ^(H) s _(i)   Formula 1

In Formula 1, F_(N) denotes an N-point FFT matrix in which each element equals $\frac{1}{\sqrt{N}}{\mathbb{e}}^{{- j}\quad 2\pi\quad{{nk}/N}}$ (n: row index, k: column index), and [•]^(H)=([•])^(T))*.

Then, L guard interval samples 20 (prefix) are attached in front of an i-th OFDM symbol block u_(i) by copy, and thus P(=L+N) number of OFDM symbol blocks u_(i,cp) 30 are formed. Here, the samples within the guard intervals satisfy the following Formula 2. u _(i,cp) (−L+n)= u _(i,cp) (N−L+n), n=0, . . . , L−1  Formula 2

However, the OFDM system has come to use a pilot for channel estimation, resulting in deterioration of frequency band efficiency. For example, among channel allocation methods of Institute of Electrical and Electronics Engineers (IEEE) 802.16e, for which standardization has recently been completed, a partial usage of subchannels (PUSC) scheme uses one pilot per six carriers to obtain desired performance, thus deteriorating system efficiency.

Since the number of pilots increases in proportion to a desired throughput, frequency band efficiency consequently undergoes severe deterioration.

Therefore, blind channel estimation algorithms are proposed as a method for increasing frequency band efficiency, the blind channel estimation algorithms estimating a channel not by using a pilot but by using an inserted cyclic prefix within a guard interval repeated in OFDM symbol blocks.

The above-described blind channel estimation methods according to conventional art have the following problems.

First, the blind channel estimation methods that have been disclosed so far are very sensitive to channel zero and have unstable characteristics. Thus, stable, high-accuracy channel estimation is not possible.

Second, since a somewhat stable and accurate blind channel estimation method would be exceedingly complex, such a method would be difficult to practically apply to blind channel estimation.

SUMMARY OF THE INVENTION

The present invention is directed to a stable channel estimation method increasing frequency band efficiency, which is deteriorated by using a pilot, and reducing complexity and sensitivity to channel zero.

The present invention is also directed to a method for directly estimating a channel from two cross matrixes obtained using a cyclic prefix of a guard interval, to increase a degree of freedom for channel zero.

The present invention is also directed to a method of stably estimating impulse response of a channel using equal gain combining and singular value decomposition (SVD) on a rank-one matrix without the influence of channel zero.

One aspect of the present invention provides a method for blind channel estimation, comprising the steps of: (a) generating an i-th symbol block S_(i) including N carriers; (b) performing an inverse fast Fourier transform (IFFT) operation on the i-th symbol block and forming an orthogonal frequency division multiplexing (OFDM) symbol block; (c) attaching a guard interval sample in front of the i-th OFDM symbol block U_(i) and forming at least one OFDM symbol block U_(i,cp); and (d) modeling the formed OFDM symbol block U_(i,cp) with a channel finite impulse response (FIR) filter h and noise v and estimating channel impulse response using signals yi received through a channel.

Step (b) may be performed by applying the following Formula 1: u _(i) =[u _(i)(0), . . . , u _(i)(N−1)]^(T) =F _(N) ^(H) s _(i)   Formula 1

In step (c), P(=L+N) number of OFDM symbol blocks U_(i,cp) may be formed, and P=qL and N=(q−1)L (q is a positive integer larger than 2).

The signals yi received in step (d) may be defined by the following Formula 4: y _(i) =H[ u ^(T) _(i−1,q−1) u ^(T) _(i,cp) ]^(T)+ v _(i)   Formula 4

The channel matrix H may be an upper or lower triangular Toeplitz matrix having a size of P×(P+L).

The length of the guard interval may be designed to be longer than the length of the channel.

The guard interval may be a vector having a length of 2L and defined by signals generated by the following Formula 5 and 6: g_(i,0) (0)=[u^(T) _(i−1,q−1) u^(T) _(i,0) ]^(T) g _(i,0) (L−1)=[u _(i−1,q−1)(L−1) u ^(T) _(i,0) u _(i,1)(0) . . . u _(i,1)(L−2)]^(T)  Formula 5 and g_(i,q−1) (0)=[u^(T) _(i,q−2) u^(T) _(i,q−1) ]^(T) g _(i,q−1) (L−1)=[u _(i,q−2)(L−1) u ^(T) _(i,q−1) u _(i+!,0)(0) . . . u _(i+1,0)(L−2)]^(T)  Formula 6

In step (d), the signals generated by Formulae 5 and 6 may pass through the channel and be defined by the following Formula 7: y _(i,0) (d)=[H ₂ H ₁] g _(i,0) (d)+[v _(i)(d) . . . v _(i)(d+L−1)]^(T) y _(i,q−1) (d)=[H ₂ H ₁] g _(i,q−1) (d)+[v _(i)((q−1)L+d) . . . v _(i)(qL−1+d)]^(T)  Formula 7 and the channel impulse response may be estimated.

The channel impulse response may be estimated by a cross-correlation matrix obtained by performing cross-correlation on the signals defined by Formula 7.

The cross-correlation matrix R₀ may be obtained by the following Formula 8: $\begin{matrix} {\begin{matrix} {R_{0} = {E\left\{ {{y_{i,0}(0)}{y_{i,{q - 1}}(0)}^{H}} \right\}}} \\ {= {H_{1}E\left\{ {u_{i,0}u_{i,{q - 1}}^{H}} \right\} H_{1}^{H}}} \\ {= {H_{1}H_{1}^{H}}} \end{matrix}\quad} & {{Formula}\quad 8} \end{matrix}$

The cross-correlation matrix R₀ may be obtained by the following Formula 9: $\begin{matrix} \begin{matrix} {R_{q - 1} = {E\left\{ {{y_{i,0}\left( {L - 1} \right)}{y_{i,{q - 1}}\left( {L - 1} \right)}^{H}} \right\}}} \\ {= {\left\lbrack {H_{2}P\quad H_{1}\underset{\_}{e_{1}}} \right\rbrack E{\left\{ {u_{i,0}u_{i,{q - 1}}^{H}} \right\}\left\lbrack {H_{2}P\quad H_{1}\underset{\_}{e_{1}}} \right\rbrack}^{H}}} \end{matrix} & {{Formula}\quad 9} \end{matrix}$

The unit column vector e₁ may have a first element that is 1 and other elements that are all 0.

Step (d) may include the steps of: (d1) calculating cross-correlations on the signals yi received through the channel; (d2) performing equal gain combining on the calculated cross-correlations and detecting at least two singular vectors by a rank-one matrix; (d3) performing singular value decomposition (SVD) on the obtained rank-one matrix and detecting the largest singular vector; and (d4) identifying a phase using a differential modulation method or one pilot carrier on the detected unit vector and estimating the channel impulse response.

In step (d1), the cross-correlations may be calculated using the following Formula 10: $\begin{matrix} {{R_{0} \approx {\frac{1}{B}{\sum\limits_{i = 0}^{B - 1}{{\underset{\_}{y_{i,0}}(0)}{\underset{\_}{y_{i,{q - 1}}}(0)}^{H}}}}}{R_{q - 1} \approx {\frac{1}{B}{\sum\limits_{i = 0}^{B - 1}{{\underset{\_}{y_{i,0}}\left( {L - 1} \right)}{\underset{\_}{y_{i,{q - 1}}}\left( {L - 1} \right)}^{H}}}}}} & {{Formula}\quad 10} \end{matrix}$

In step (d2), the rank-one matrix H may be obtained using the following Formula 19: $\begin{matrix} {{\hat{h}{\hat{h}}^{H}} = \frac{R_{0,h} + R_{{q - 1},h}}{2}} & {{Formula}\quad 19} \end{matrix}$

Step (d2) may include the steps of: directly calculating a first estimated singular vector from a first column of a matrix R₀ of the calculated cross-correlations; and directly calculating a second estimated singular vector from a last row of a matrix R_(q−1) of the calculated cross-correlations.

The first estimated unit vector may be calculated by the following Formulae 11 and 12: $\begin{matrix} {{{1{st}\quad{column}\quad{of}\quad R_{0}} = {\underset{\_}{r_{1}} = \left( {\begin{bmatrix} h_{0} & h_{1} & \ldots & h_{L - 1} \end{bmatrix}h_{0}^{*}} \right)^{T}}}{and}} & {{Formula}\quad 11} \\ {{\hat{h}}_{1} = \frac{\underset{\_}{r_{1}}}{\sqrt{{first}\quad{element}\quad{of}\quad\underset{\_}{r_{1}}}}} & {{Formula}\quad 12} \end{matrix}$

The second estimated unit vector may be calculated by the following Formulae 13 and 14: $\begin{matrix} {{{{last}\quad{row}\quad{of}\quad R_{q - 1}} = {\overset{\_}{r_{L}} = {h_{L - 1}\begin{bmatrix} h_{0}^{*} & h_{1}^{*} & \ldots & h_{L - 1}^{*} \end{bmatrix}}}}{and}} & {{Formula}\quad 13} \\ {{\hat{h}}_{2} = \left( \frac{\overset{\_}{r_{L}}}{\sqrt{{first}\quad{element}\quad{of}\quad\overset{\_}{r_{L}}}} \right)^{H}} & {{Formula}\quad 14} \end{matrix}$

When h₀, h_(q−1) satisfies h₀,h_(q−1)≠0 and h₀,h_(q−1)>>0, the equal gain combining may be performed by the following Formula 15: $\begin{matrix} {\hat{h} = {\frac{{\hat{h}}_{1} + {\hat{h}}_{2}}{2} = {\left( \frac{{h_{0}^{*}/{h_{0}}} + {h_{L - 1}^{*}/{h_{L - 1}}}}{2} \right)\begin{bmatrix} h_{0} & h_{1} & \ldots & h_{L - 1} \end{bmatrix}}}} & {{Formula}\quad 15} \end{matrix}$

In step (d3), the largest singular vector may be detected by at least two iterations of an SVD algorithm.

Step (d) may include the steps of: calculating cross-correlations on the signals yi received through the channel; detecting a singular vector by a lower Toeplitz matrix obtained by performing Cholesky factorization on the calculated cross-correlations; and identifying a phase using a differential modulation method or one pilot carrier on the detected unit vector and estimating the channel impulse response.

BRIEF DESCRIPTION OF THE DRAWINGS

The above and other features and advantages of the present invention will become more apparent to those of ordinary skill in the art by describing in detail exemplary embodiments thereof with reference to the attached drawings in which:

FIG. 1 is a diagram illustrating a method for estimating a channel by inserting a cyclic prefix into a guard interval in an orthogonal frequency division multiplexing (OFDM) system;

FIG. 2 is a diagram illustrating a method for estimating a channel by inserting a cyclic prefix into a guard interval in an OFDM system according to an exemplary embodiment of the present invention;

FIG. 3 illustrates a Toeplitz structure of a channel matrix within a guard interval in an OFDM system according to an exemplary embodiment of the present invention;

FIG. 4 illustrates the structure of a lower triangular Toeplitz matrix within a guard interval in an OFDM system according to an exemplary embodiment of the present invention;

FIG. 5 illustrates the structure of a upper triangular Toeplitz matrix within a guard interval in an OFDM system according to an exemplary embodiment of the present invention; and

FIG. 6 illustrates the structure of a left shifting matrix within a guard interval in an OFDM system according to an exemplary embodiment of the present invention.

DETAILED DESCRIPTION OF EXEMPLARY EMBODIMENTS

Hereinafter, exemplary embodiments of the present invention will be described in detail. However, the present invention is not limited to the exemplary embodiments disclosed below and can be implemented in various forms. Therefore, the present exemplary embodiments are provided for complete disclosure of the present invention and to fully convey the scope of the present invention to those of ordinary skill in the art.

A method for blind channel estimation using a guard interval of a signal in an orthogonal frequency division multiplexing (OFDM) system according to an exemplary embodiment of the present invention will be described with reference to the appended drawings. Here, as illustrated in FIG. 1, the length of a guard interval in an actual OFDM system is generally defined as the total fast Fourier transform (FFT) length N divided by an integer (N/a, a is an integer). Therefore, for ease of description, the present invention assumes that P=qL and N=(q−1)L (q is a positive integer larger than 2).

In this connection, u_(i) 30 of FIG. 1 can be divided into a total number q of sub-blocks each having a length of L, as defined by the following Formula 3, which is illustrated in FIG. 2. u_(i,cp) ^(T) =[u_(i,0) ^(T) . . . u_(i,q−1) ^(T)]^(T) u _(i,b) ^(T) =[u _(i)(bL) . . . u _(i)(bL+L−1)]^(T) , b=0, . . . , q−1  Formula 3

Since a guard interval is generated by cyclic prefixing in Formula 3, a sub-block u_(i,0) =u_(i,q−1) 500.

In addition, a channel is modeled with a channel finite impulse response (FIR) filter h=[h₀ . . . h_(L−1)]^(T) and noise v.

Here, it is assumed that the length of the channel FIR is M, and the guard interval of the OFDM system is designed to be longer than the channel length, i.e., M≦L. This means h₁=0, 1≧M.

Therefore, the channel is a Toeplitz matrix having a size of P×(P+L) illustrated in FIG. 3, in which a first column corresponds to [0 . . . 0]^(T), and a first row corresponds to [0 h_(L−1) . . . h₀ 0 . . . 0].

Consequently, an i-th block received through the channel may be expressed by the following Formula 4. y _(i) =H[ u ^(T) _(i−1,q−1) u ^(T) _(i,cp) ]^(T)+ v _(i)   Formula 4

In Formula 4, u_(i−1,q−1) denotes a signal of a (q−1)th sub-block of an (i−1)th transmission block, v_(i) =[v_(i)(0) . . . v_(i)(qL−1)]^(T) denotes an i-th block noise signal vector, and a sample v_(i)(n) denotes complex Gaussian white noise having a zero average and variance of noise power σ².

The method for blind channel estimation proposed as described above fundamentally utilizes information on a guard interval repeated in an OFDM symbol block. To this end, a lower triangular Toeplitz matrix H₁ and an upper triangular Toeplitz matrix H₂ illustrated in FIGS. 4 and 5 are defined.

In addition, a vector having a guard interval length of 2L is defined by signals 100 and 200 generated by the following Formula 5 to transmit the vector. g_(i,0) (0)=[u^(T) _(i−1,q−1) u^(T) _(i,0) ]^(T) g _(i,0) (L−1)=[u _(i−1,q−1)(L−1) u ^(T) _(i,0) u _(i,1)(0) . . . u _(i,1)(L−2)]^(T)  Formula 5

Likewise, signals 300 and 400 are generated by the following Formula 6. g_(i,q−1) (0)=[u^(T) _(i,q−2) u^(T) _(i,q−1) ]^(T) g _(i,q−1) (L−1)=[u _(i,q−2)(L−1) u ^(T) _(i,q−1) u _(i+!,0)(0) . . . u _(i+1,0)(L−2)]^(T)  Formula 6

The signals of Formulae 5 and 6 each pass through the channel and may be expressed by the following Formula 7 in case of d=0, L−1. y _(i,0) (d)=[H ₂ H ₁] g _(i,0) (d)+[v _(i)(d) . . . v _(i)(d+L−1)]^(T) y _(i,q−1) (d)=[H ₂ H ₁] g _(i,q−1) (d)+[v _(i)((q−1)L+d) . . . v _(i)(qL−1+d)]^(T)  Formula 7

The method for blind channel estimation is based on generating a correlation matrix by performing cross-correlation on y_(i,0) (d) and y_(i,q−1) (d). Assuming that a cross-correlation matrix is R₀ when d=0, and R_(q−1) when d=q−1, expected values of the cross-correlation matrix obtained as a result of performing the cross-correlation may be expressed by the following Formulae 8 and 9. $\begin{matrix} {\begin{matrix} {R_{0} = {E\left\{ {{y_{i,0}(0)}{y_{i,{q - 1}}(0)}^{H}} \right\}}} \\ {= {H_{1}E\left\{ {u_{i,0}u_{i,{q - 1}}^{H}} \right\} H_{1}^{H}}} \\ {= {H_{1}H_{1}^{H}}} \end{matrix}\quad} & {{Formula}\quad 8} \\ \begin{matrix} {R_{q - 1} = {E\left\{ {{y_{i,0}\left( {L - 1} \right)}{y_{i,{q - 1}}\left( {L - 1} \right)}^{H}} \right\}}} \\ \left. {\left. {= \begin{matrix} \left\lbrack {H_{2}P} \right. & {H_{1}\underset{\_}{e_{1}}} \end{matrix}} \right\rbrack E\left\{ {u_{i,0}u_{i,{q - 1}}^{H}} \right\}\begin{matrix} \left\lbrack {H_{2}P} \right. & {H_{1}\underset{\_}{e_{1}}} \end{matrix}} \right\rbrack^{H} \end{matrix} & {{Formula}\quad 9} \end{matrix}$

Here, P is a left shifting matrix having a size of L×L and is defined as illustrated in FIG. 6. In addition, e₁ is a unit column vector having a length of L, a first element that is 1, and other elements that are all 0.

As shown in Formulae 8 and 9, when a correlation matrix is generated by performing cross-correlation rather than auto-correlation, it can be seen that influence of noise power σ² is theoretically eliminated. Thus, the method becomes strong against noise due to the cross-correlation.

In practical application, R₀, R_(q−1) can be estimated by averaging B symbol blocks as in the following Formula 10. $\begin{matrix} {{R_{0} \approx {\frac{1}{B}{\sum\limits_{i = 0}^{B - 1}{{\underset{\_}{y_{i,0}}(0)}{\underset{\_}{y_{i,{q - 1}}}(0)}^{H}}}}}{R_{q - 1} \approx {\frac{1}{B}{\sum\limits_{i = 0}^{B - 1}{{\underset{\_}{y_{i,0}}\left( {L - 1} \right)}{\underset{\_}{y_{i,{q - 1}}}\left( {L - 1} \right)}^{H}}}}}} & {{Formula}\quad 10} \end{matrix}$

A variety of methods are available for estimating impulse response of the channel from R₀, R_(q−1) given by Formula 10.

An estimated channel impulse value obtained by the following methods in the end is defined as ĥ.

The first method most simply estimates channel impulse response from given R₀, R_(q−1). This includes a method of directly obtaining an estimation value ĥ₁ from the first column of R₀ and a method of directly obtaining an estimation value ĥ₂ from the last row of R_(q−1). Since the first column of R₀ is defined by Formula 11, the method can estimate ĥ₁ by the following Formula 12 without specifying phase. $\begin{matrix} {{1\quad{st}\quad{column}\quad{of}\quad R_{0}} = {\underset{\_}{r_{1}} = \left( {\begin{bmatrix} h_{0} & h_{1} & \ldots & h_{L - 1} \end{bmatrix}h_{0}^{*}} \right)^{T}}} & {{Formula}\quad 11} \\ {{\hat{h}}_{1} = \frac{\underset{\_}{r_{1}}}{\sqrt{{first}\quad{element}\quad{of}\quad\underset{\_}{r_{1}}}}} & {{Formula}\quad 12} \end{matrix}$

Likewise, since the last row of R_(q−1) is defined by Formula 13, the method can estimate ĥ₂ by the following Formula 14 without specifying phase. $\begin{matrix} {{{last}\quad{row}\quad{of}\quad R_{q - 1}} = {\overset{\_}{r_{L}} = {h_{L - 1}\begin{bmatrix} h_{0}^{*} & h_{1}^{*} & \ldots & h_{L - 1}^{*} \end{bmatrix}}}} & {{Formula}\quad 13} \\ {{\hat{h}}_{2} = \left( \frac{\overset{\_}{r_{L}}}{\sqrt{{last}\quad{element}\quad{of}\quad\overset{\_}{r_{L}}}} \right)^{H}} & {{Formula}\quad 14} \end{matrix}$

The vagueness of a phase existing in Formulae 12 and 14 can be eliminated using a differential modulation method or one pilot carrier at a transmitting end.

Here, when the result of Formula 12 is combined with that of Formula 14, it is possible to further increase the accuracy of estimated channel impulse response.

Then, equal gain combining is performed on the results of Formulae 12 and 14, and a result corresponding to following Formula 15 can be obtained. $\begin{matrix} {\hat{h} = {\frac{{\hat{h}}_{1} + {\hat{h}}_{2}}{2} = {\left( \frac{\frac{h_{0}^{*}}{h_{0}} + \frac{h_{L - 1}^{*}}{h_{L - 1}}}{2} \right)\begin{bmatrix} h_{0} & h_{1} & \ldots & h_{L - 1} \end{bmatrix}}}} & {{Formula}\quad 15} \end{matrix}$

However, the first method using Formulae 11 to 15 must satisfy the assumption that h₀, h_(q−1) are not 0. In addition, when one of h₀, h_(q−1) is almost 0, an error is greatly amplified while Formulae 12 and 14 are calculated. Consequently, inaccurate channel impulse response is estimated by Formula 15.

Therefore, the first method can be used when the following Formula 16 is satisfied. h₀, h_(q−1)≠0 and h₀, h_(q−1)>>0  Formula 16

When only h₀ satisfies Formula 16, the method according to Formulae 11 and 12 must be used, and when only h_(q−1) satisfies Formula 16, the method according to Formulae 13 and 14 must be used.

Another available method obtains ĥ using Cholesky factorization.

Since H₁ is a lower triangular Toeplitz matrix as shown in FIG. 4, it is possible to easily estimate channel impulse response by Cholesky factorization of R₀.

More specifically, when the Cholesky factorization of R₀, which is approximately obtained by Formula 10, is performed, a lower triangular Toeplitz matrix {tilde over (H)}₁ (estimated values for H₁) approximate to H₁ can be obtained. In this case, the first column of {tilde over (H)}₁ is reflected in ĥ with unspecified phase. However, in order to apply Cholesky factorization, R₀ must be always positive definite.

Even when little noise exists in R₀, R₀ becomes non-positive definite, and thus Cholesky factorization diverges. Therefore, the Cholesky method cannot be used in actual systems.

The second method stably obtains an estimated channel value ĥ using R₀, R_(q−1) at the same time. This requires neither the condition of Formula 16 nor the condition that R₀ must be positive definite.

The second method proposed by the present invention obtains ĥ by the following Formulae 17 and 18. R _(0,h) =R ₀ −PR ₀ P ^(T)  Formula 17 R _(q−1,h) =R _(q−1) −P ^(T) R _(q−1) P  Formula 18

Then, it is possible to generate a matrix corresponding to the outer product of ĥ and ĥ^(H) by performing equal gain combining on results obtained by Formulae 17 and 18 as shown in the following Formula 19. In other words, a rank-one matrix can be generated. $\begin{matrix} {{\hat{h}{\hat{h}}^{H}} = \frac{R_{0,h} + R_{{q - 1},h}}{2}} & {{Formula}\quad 19} \end{matrix}$

In order to describe this method in detail, it is assumed that a j-th column of H₁ is h_(j) ¹ as shown in FIG. 4 and that a j-th column of H₂ is h_(j) ^(q−1) as shown in FIG. 5.

Then, Formula 17 can be combined with Formula 8 to yield the following Formula 20. $\begin{matrix} \begin{matrix} {R_{0,h} = {{\begin{bmatrix} \underset{\_}{h_{0}^{0}} & \underset{\_}{h_{1}^{0}} & \underset{\_}{h_{2}^{0}} & \ldots & \underset{\_}{h_{L - 1}^{0}} \end{bmatrix}\begin{bmatrix} \underset{\_}{h_{0}^{0}} & \underset{\_}{h_{1}^{0}} & \underset{\_}{h_{2}^{0}} & \ldots & \underset{\_}{h_{L - 1}^{0}} \end{bmatrix}}^{H} -}} \\ {{{P\begin{bmatrix} \underset{\_}{h_{0}^{0}} & \underset{\_}{h_{1}^{0}} & \underset{\_}{h_{2}^{0}} & \ldots & \underset{\_}{h_{L - 1}^{0}} \end{bmatrix}}\begin{bmatrix} \underset{\_}{h_{0}^{0}} & \underset{\_}{h_{1}^{0}} & \underset{\_}{h_{2}^{0}} & \ldots & \underset{\_}{h_{L - 1}^{0}} \end{bmatrix}}^{H}P^{T}} \\ {= {{\underset{\_}{h_{0}^{0}}\quad\underset{\_}{h_{0}^{0^{H}}}} + {\underset{\_}{h_{1}^{0}}\quad\underset{\_}{h_{1}^{0^{H}}}} + \ldots + {\underset{\_}{h_{L - 1}^{0}}\quad\underset{\_}{h_{L - 1}^{0^{H}}}} - {P\quad\underset{\_}{h_{0}^{0}}\quad\underset{\_}{h_{0}^{0^{H}}}P^{T}} -}} \\ {{P\quad\underset{\_}{h_{1}^{0}}\quad\underset{\_}{h_{1}^{0^{H}}}P^{T}} - \ldots - {P\quad\underset{\_}{h_{L - 1}^{0}}\quad\underset{\_}{h_{L - 1}^{0^{H}}}P^{T}}} \end{matrix} & {{Formula}\quad 20} \end{matrix}$

The following Formula 21 can be deduced from Formula 20. $\begin{matrix} \begin{matrix} {P\underset{\_}{h_{0}^{0}}\quad\underset{\_}{h_{0}^{0^{H}}}P^{T}} & = & {\underset{\_}{h_{1}^{0}}\quad\underset{\_}{h_{1}^{0^{H}}}} \\ {P\underset{\_}{h_{1}^{0}}\quad\underset{\_}{h_{1}^{0^{H}}}P^{T}} & = & {\underset{\_}{h_{2}^{0}}\quad\underset{\_}{h_{2}^{0^{H}}}} \\ \vdots & = & \vdots \\ {P\underset{\_}{h_{L - 2}^{0}}\quad\underset{\_}{h_{L - 2}^{0^{H}}}P^{T}} & = & {\underset{\_}{h_{L - 1}^{0}}\quad\underset{\_}{h_{L - 1}^{0^{H}}}} \\ {P\underset{\_}{h_{L - 1}^{0}}\quad\underset{\_}{h_{L - 1}^{0^{H}}}P^{T}} & = & 0 \end{matrix} & {{Formula}\quad 21} \end{matrix}$

Consequently, by applying Formula 21, only h₀ ⁰ h₀ ⁰ ^(H) remains in Formula 20.

In the same manner, Formula 18 can be combined with Formula 9 to yield the following Formula 22. $\begin{matrix} \begin{matrix} {R_{{q - 1},h} = \begin{bmatrix} h_{1}^{q - 1} & h_{2}^{q - 1} & h_{3}^{q - 1} & \ldots & h_{L - 1}^{q - 1} & h_{0}^{0} \end{bmatrix}} \\ {\begin{bmatrix} h_{1}^{q - 1} & h_{2}^{q - 1} & h_{3}^{q - 1} & \ldots & h_{L - 1}^{q - 1} & h_{0}^{0} \end{bmatrix}^{H} -} \\ {P^{T}\begin{bmatrix} h_{1}^{q - 1} & h_{2}^{q - 1} & h_{3}^{q - 1} & \ldots & h_{L - 1}^{q - 1} & h_{0}^{0} \end{bmatrix}} \\ {\begin{bmatrix} h_{1}^{q - 1} & h_{2}^{q - 1} & h_{3}^{q - 1} & \ldots & h_{L - 1}^{q - 1} & h_{0}^{0} \end{bmatrix}^{H}P} \\ {= {{h_{1}^{q - 1}\quad h_{1}^{q - 1^{H}}} + {h_{2}^{q - 1}\quad h_{2}^{q - 1^{H}}} + \ldots +}} \\ {{h_{L - 1}^{q - 1}\quad h_{L - 1}^{q - 1^{H}}} + {h_{0}^{0}\quad h_{0}^{0^{H}}} - {{P\quad}^{T}h_{1}^{q - 1}\quad h_{1}^{q - 1^{H}}P} -} \\ {{{P\quad}^{T}h_{2}^{q - 1}\quad h_{2}^{q - 1^{H}}P} - \ldots -} \\ {{P^{T}\quad h_{L - 1}^{q - 1}\quad h_{L - 1}^{q - 1^{H}}P} - {P^{T}h_{0}^{0}h_{0}^{0^{H}}P}} \end{matrix} & {{Formula}\quad 22} \end{matrix}$

The following Formula 23 can be deduced from Formula 22. $\begin{matrix} \begin{matrix} {P^{T}h_{1}^{q - 1}\quad h_{1}^{q - 1^{H}}P} & = & 0 \\ {P^{T}h_{2}^{q - 1}\quad h_{2}^{q - 1^{H}}P} & = & {h_{1}^{q - 1}\quad h_{1}^{q - 1^{H}}} \\ \vdots & \quad & \vdots \\ {P^{T}h_{L - 1}^{q - 1}\quad h_{L - 1}^{q - 1^{H}}P} & = & {h_{L - 2}^{q - 1}\quad h_{L - 2}^{q - 1^{H}}} \\ {P^{T}h_{0}^{0}\quad h_{0}^{0^{H}}P} & = & {h_{L - 1}^{q - 1}\quad h_{L - 1}^{q - 1^{H}}} \end{matrix} & {{Formula}\quad 23} \end{matrix}$

Consequently, by applying Formula 23, only h₀ ⁰ h₀ ⁰ ^(H) remains in Formula 22.

When equal gain combining is performed on the results obtained by Formulae 17 and 18 as defined in Formula 19, it is possible to estimate h₀ ₀ h₀ ⁰ ^(H) and also reduce a caused error.

Here, in order to obtain ĥ from the result (rank-one matrix) of Formula 19, singular value decomposition (SVD) is performed.

More specifically, the matrix of Formula 19 is a square matrix and a Hermitian matrix at the same time, and thus has one main (large) singular value and other values of 0. Therefore, when an iterative SVD algorithm, e.g., Jocobi SVD, is used, a singular vector corresponding to the largest singular value can be obtained by only 1 or 2 iterations.

The obtained singular vector becomes ĥ.

The vagueness of a phase in ĥ obtained in this manner can be eliminated using the differential modulation method or one pilot carrier at a transmitting end.

The first of the above-described two methods estimates channel impulse response by Formulae 11 to 15 when Formula 16 is satisfied. In addition, the method according to Formulae 11 and 12 must be used when only h₀ satisfies Formula 16, and the method according to Formulae 13 and 14 must be used when only h_(q−1) satisfies Formula 16.

On the contrary, when any of the above-mentioned conditions is not satisfied, the first method cannot be used.

In addition, the second method can reduce complexity and stably estimate a channel with high accuracy.

As described above, a method for blind channel estimation using a guard interval of a signal in an OFDM technique according to the present invention does not require any change in conventional OFDM transmitters. The proposed methods are low-complexity blind channel estimation algorithms that can be applied to actual systems, do not affect channel zero, and are stable. Also, the methods can be modified into semi-blind methods.

In addition, the SVD of a rank-one matrix has low complexity in practical application and thus can be used for actual blind channel estimation.

While the invention has been shown and described with reference to certain exemplary embodiments thereof, it will be understood by those skilled in the art that various changes in form and details may be made therein without departing from the spirit and scope of the invention as defined by the appended claims. 

1. A method for blind channel estimation, comprising the steps of: (a) generating an i-th symbol block S_(i) including N carriers; (b) performing an inverse fast Fourier transform (IFFT) operation on the i-th symbol block and forming an orthogonal frequency division multiplexing (OFDM) symbol block; (c) attaching a guard interval sample in front of the i-th OFDM symbol block U_(i) and forming at least one OFDM symbol block U_(i,cp); and (d) modeling the formed OFDM symbol block U_(i,cp) with a channel finite impulse response (FIR) filter h and noise v and estimating channel impulse response using signals yi received through a channel, wherein N, i, L are non-zero positive integers.
 2. The method of claim 1, wherein step (b) is performed by the following Formula 1: u _(i) =[u _(i)(0), . . . , u _(i)(N−1)]^(T) =F _(N) ^(H) s _(i) ,  Formula 1wherein F_(N) is an N-point fast Fourier transform (FFT) matrix.
 3. The method of claim 1, wherein in step (c), P(=L+N) number of OFDM symbol blocks U_(i,cp) are formed, and P=qL and N=(q−1)L, wherein q is a positive integer larger than
 2. 4. The method of claim 1, wherein the signals yi received in step (d) are defined by the following Formula 4: y _(i) =H[ u ^(T) _(i−1,q−1) u ^(T) _(i,cp) ]^(T)+ v _(i) ,  Formula 4wherein H is a channel matrix, u_(i−1,q−1) is a signal of a (q−1)th sub-block of an (i−1)th transmission block, and v_(i) is an i-th block noise signal vector.
 5. The method of claim 4, wherein the channel matrix H is an upper or lower triangular Toeplitz matrix having a size of P×(P+L).
 6. The method of claim 1, wherein the length of the guard interval is designed to be longer than the length of the channel.
 7. The method of claim 1, wherein the guard interval is a vector having a length of 2L and defined by signals generated by the following Formulae 5 and 6: g _(i,0) (0)=[ u ^(T) _(i−1,q−1) u ^(T) _(i,0) ]^(T) g _(i,0) (L−1)=[u _(i−1,q−1)(L−1) u ^(T) _(i,0) u _(i,1)(0) . . . u _(i,1)(L−2)]^(T)  Formula 5 g _(i,q−1) (0)=[ u ^(T) _(i,q−2) u ^(T) _(i,q−1) ]^(T) g _(i,q−1) (L−1)=[u _(i,q−2)(L−1) u ^(T) _(i,q−1) u _(i+!,0)(0) . . . u _(i+1,0)(L−2)]^(T)  Formula 6
 8. The method of claim 7, wherein in step (d), the signals generated by Formulae 5 and 6 pass through the channel and are defined by the following Formula 7: y _(i,0) (d)=[H ₂ H ₁] g _(i,0) (d)+[v _(i)(d) . . . v _(i)(d+L−1)]^(T) y _(i,q−1) (d)=[H ₂ H ₁] g _(i,q−1) (d)+[v _(i)((q−1)L+d) . . . v _(i)(qL−1+d)]^(T)  Formula 7 and the channel impulse response is estimated.
 9. The method of claim 8, wherein the estimation of the channel impulse response is made by a cross-correlation matrix obtained by performing cross-correlation on the signals defined by Formula
 7. 10. The method of claim 9, wherein the cross-correlation matrix R₀ is obtained using the following Formula 8: $\begin{matrix} \begin{matrix} {R_{0} = {E\left\{ {{y_{i,0}(0)}{y_{i,{q - 1}}(0)}^{H}} \right\}}} \\ {= {H_{1}E\left\{ {u_{i,0}u_{i,{q - 1}}^{H}} \right\} H_{1}^{H}}} \\ {= {H_{1}{H_{1}^{H}.}}} \end{matrix} & {{Formula}\quad 8} \end{matrix}$
 11. The method of claim 9, wherein the cross-correlation matrix R₀ is obtained using the following Formula 9: $\begin{matrix} \begin{matrix} {R_{q - 1} = {E\left\{ {{y_{i,0}\left( {L - 1} \right)}{y_{i,{q - 1}}\left( {L - 1} \right)}^{H}} \right\}}} \\ {{= {\left\lbrack {H_{2}P\quad H_{1}\underset{\_}{e_{1}}} \right\rbrack E{\left\{ {u_{i,0}u_{i,{q - 1}}^{H}} \right\}\left\lbrack {H_{2}P\quad H_{1}\underset{\_}{e_{1}}} \right\rbrack}^{H}}},} \end{matrix} & {{Formula}\quad 9} \end{matrix}$ wherein P is a left shifting matrix having a size of L×L, and e₁ is a unit column vector having a length of L.
 12. The method of claim 11, wherein the unit column vector has a first element that is 1 and other elements that are all
 0. 13. The method of claim 1, wherein step (d) comprises the steps of: (d1) calculating cross-correlations on the signals yi received through the channel; (d2) performing equal gain combining on the calculated cross-correlations and detecting at least two singular vectors by a rank-one matrix; (d3) performing a singular value decomposition (SVD) on the obtained rank-one matrix and detecting a largest singular vector; and (d4) identifying a phase using a differential modulation method or one pilot carrier on the detected unit vector and estimating the channel impulse response.
 14. The method of claim 13, wherein in step (d1), the cross-correlations are calculated using the following Formula 10: $\begin{matrix} {{R_{0} \approx {\frac{1}{B}{\sum\limits_{i = 0}^{B - 1}{{\underset{\_}{y_{i,0}}(0)}{\underset{\_}{y_{i,{q - 1}}}(0)}^{H}}}}}{R_{q - 1} \approx {\frac{1}{B}{\sum\limits_{i = 0}^{B - 1}{{\underset{\_}{y_{i,0}}\left( {L - 1} \right)}{\underset{\_}{y_{i,{q - 1}}}\left( {L - 1} \right)}^{H}}}}}} & {{Formula}\quad 10} \end{matrix}$
 15. The method of claim 13, wherein in step (d2), the rank-one matrix H is obtained using the following Formula 19: $\begin{matrix} {{{\hat{h}\quad{\hat{h}}^{H}} = \frac{R_{0,h} + R_{{q - 1},h}}{2}},} & {{Formula}\quad 19} \end{matrix}$ wherein R_(0,h)=R₀−PR₀P^(T) and R_(q−1,h)=R_(q−1)−P^(T)R_(q−1)P.
 16. The method of claim 13, wherein step (d2) comprises the steps of: directly calculating a first estimated singular vector from a first column of a matrix R₀ of the calculated cross-correlations; and directly calculating a second estimated singular vector from a last row of a matrix R_(q−1) of the calculated cross-correlations.
 17. The method of claim 16, wherein the first estimated unit vector is calculated by the following Formulae 11 and 12: $\begin{matrix} {{1{st}\quad{column}\quad{of}\quad R_{0}} = {\underset{\_}{r_{1}} = {\left( {\begin{bmatrix} h_{0} & h_{1} & \cdots & h_{L - 1} \end{bmatrix}h_{0}^{*}} \right)^{T}{and}}}} & {{Formula}\quad 11} \\ {{\hat{h}}_{1} = \frac{\underset{\_}{r_{1}}}{\sqrt{{first}\quad{element}\quad{of}\quad\underset{\_}{r_{1}}}}} & {{Formula}\quad 12} \end{matrix}$
 18. The method of claim 16, wherein the second estimated unit vector is calculated by the following Formulae 13 and 14: $\begin{matrix} {{{last}\quad{row}\quad{of}\quad R_{q - 1}} = {\overset{\_}{r_{L}} = {{h_{L - 1}\begin{bmatrix} h_{0}^{*} & h_{1}^{*} & \cdots & h_{L - 1}^{*} \end{bmatrix}}\quad{and}}}} & {{Formula}\quad 13} \\ {h_{2} = \left( \frac{\overset{\_}{r_{L}}}{\sqrt{{last}\quad{element}\quad{of}\quad\overset{\_}{r_{L}}}} \right)^{H}} & {{Formula}\quad 14} \end{matrix}$
 19. The method of claim 17 or 18, wherein when h₀, h_(q−1) satisfies h₀, h_(q−1)≠0 and h₀, h_(q−1)>>0, the equal gain combining is performed by the following Formula 15: $\begin{matrix} \begin{matrix} {\hat{h} = \frac{{\hat{h}}_{1} + {\hat{h}}_{2}}{2}} \\ {= {\left( \frac{{h_{0}^{*}/{h_{0}}} + {h_{L - 1}^{*}/{h_{L - 1}}}}{2} \right)\begin{bmatrix} h_{0} & h_{1} & \cdots & h_{L - 1} \end{bmatrix}}} \end{matrix} & {{Formula}\quad 15} \end{matrix}$
 20. The method of claim 13, wherein in step (d3), the largest singular vector is detected by at least two iterations of an SVD algorithm.
 21. The method of claim 1, wherein step (d) comprises the steps of: calculating cross-correlations on the signals yi received through the channel; detecting a singular vector by a lower Toeplitz matrix obtained by performing Cholesky factorization on the calculated cross-correlations; and identifying a phase using a differential modulation method or one pilot carrier on the detected unit vector and estimating the channel impulse response. 